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Abstract — Based on the Ogorodnikov-Milne model, we analyze the proper motions of Tycho-2 and UCAC2 stars. 
We have established that the model component that describes the rotation of all stars under consideration around the 
Galactic y axis differs significantly from zero at various magnitudes. We interpret this rotation found using the most 
distant stars as a residual rotation of the ICRS /Tycho-2 system relative to the inertial reference frame. For the most 
distant (d « 900 pc) Tycho-2 and UCAC2 stars, the mean rotation around the Galactic y axis has been found to be 
M13 = — 0.37 ± 0.04 mas yr _1 . The proper motions of UCAC2 stars with magnitudes in the range 12 — 15 m are shown 
to be distorted appreciably by the magnitude equation in /i a cos<5, which has the strongest effect for northern-sky stars 
with a coefficient of —0.60 ± 0.05 mas yr" 1 mag" 1 . We have detected no significant effect of the magnitude equation in 
the proper motions of UCAC2 stars brighter than « ll m . 
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INTRODUCTION 

The current standard International Celestial Reference 
System (ICRS) is realized by a catalog of radio posi- 
tions for 212 compact extragalactic sources uniform- 
ly distributed over the entire sky and observed by 
very long baseline interferometry (Ma et al. 1998). In 
the optical range, the first ICRS realization was the 
Hipparcos catalog (ESA 1997). Application of various 
methods of analysis reveals a small residual rotation 
of ICRS/Hipparcos relative to the inertial reference 
system with « —0.4 ± 0.1 mas yr" 1 (Bobylev 2004a, 
2004b). 

The method considered here is based on the study 
of two solid-body rotation tensor components that de- 
scribe the rotation around the y and x axes in the 
Galactic coordinate system. Application of this method 
to the proper motions of TRC stars (H0g et al. 1998) 
showed that the determination of the rotation around 
the y axis is noticeably affected by the actual rotation 
of the Local system's stars (Bobylev 2004a). In this pa- 
per, we consider the Tycho-2 catalog (H0g et al. 2000), 
which does not differ systematically from the TRC cat- 
alog, but contains twice as many stars. We expect a 
confirmation of the results obtained using TRC stars. 

Invoking fairly distant stars that are free from the 
effects of both streams close to the Sun and the Local 



system of stars as a whole is topical for a reliable appli- 
cation of the method. The UCAC2 catalog (Zacharias 
et al. 2004), which contains positions and proper mo- 
tions for some 48 million stars, is of great interest in 
this connection. This catalog extends the ICRS/Tycho- 
2 system to the 17 th magnitude. 

The goal of this paper is to study the kinematic 
parameters of a large number of stars as a function of 
their distance. To estimate the distances, we use the 
method of comparing the statistical parallaxes with the 
solar velocity (Oiling and Dehnen 2003); the latter is 
currently known well (Dehnen and Binney 1998). 

1 THE MODEL 

In this paper, we use a rectangular Galactic coordinate 
system with the axes directed away from the observer 
toward the Galactic center (I = 0°,b — 0°, the x axis 
or axis 1), along the Galactic rotation (I = 90°, 6 = 0°, 
the y axis or axis 2), and toward the North Galactic 
pole (b — 90° , the z axis or axis 3) . In the Ogorodnikov- 
Milne model, we use the notation introduced by Clube 
(1972, 1973) and employed by du Mont (1977, 1978). 
When only the stellar proper motions are used, one 
of the diagonal terms of the local deformation tensor 
is known (Ogorodnikov 1965) to remain indeterminate. 



Therefore, we determine differences of the form (M+ — 
M+) and (M3+3-M+). 

The conditional equations can be written as 

/itz cos 6 = (l/r)(X Q sinZ — Y" Q cosZ) — (1) 
— M~ 2 cos Z sin 6 — M~ sin Z sin 6 + M~ cos 6+ 
+M+ cos 2 Z cos 6 - M+ sin Z sin 6+ 
+M 23 cos Z sin 6— 
-0.5(M+ - M+) sin 2Z cos 6, 
m = (l/r) (Xq cos Z sin b + Yq sin I sin & — (2) 
— Z Q cos 6) + M~ 2 sin Z — M~ cos Z— 
-0.5M+ sin 2Z sin 26 + M+ cos Z cos 26+ 
+M+ sin Z cos 26- 
-0.5(M+ - M+) cos 2 Z sin 26+ 
+0.5(M3 + 3-M 2 + 2 )sin26. 

where Xq, Yq, Zq are the peculiar solar velocity com- 
ponents and M~ 2 , M~ 3 , M~ 3 are the components of the 
vector of solid-body rotation of a small solar neigh- 
borhood around the corresponding axes. In accordance 
with the adopted rectangular coordinate system, the 
following rotations are positive: from axis 1 to axis 2, 
from axis 2 to axis 3, and from axis 3 to axis 1. The 
quantity M 21 (mas yr _1 ) is related to the Oort con- 
stant B (km s _1 kpc -1 ) via the proportionality coef- 
ficient 4.74. Each of the quantities M+, M+ and M+ 
describes the deformation in the corresponding plane. 
The quantity M+ (mas yr _1 ) is related to the Oort 
constant A (km s _1 kpc -1 ) via the proportionality co- 
efficient 4.74. The diagonal components of the local de- 
formation tensor M+ , M+ and M 3 + describe the overall 
contraction or expansion of the entire stellar system. 
The system of conditional equations (1) and (2) in- 
cludes eleven sought-for unknowns to be determined 
by the least squares method. As can be seen from Eq. 
(1), the two pairs of unknowns M~ 3 and M+ as well as 
M~ 2 and M+ have the same coefficients, sin Z sin 6 and 
cos Z sin 6, respectively. As a result, the variables are 
ill-separated. In the UCAC2 section, we analyze the 
results of both the simultaneous solution of Eqs. (1) 
and (2) and the separate solution of only Eq. (2). The 
quantity l/r is a parallactic factor that is taken to be 
unity when solving Eqs. (1) and (2). In this case, the 
stars are referred to a unit sphere. In this approach, all 
of the parameters being determined are proportional 
to the heliocentric distance of the stellar centroid un- 
der consideration and are expressed in the same units 
as the stellar proper motion components, i.e., in mas 
yr _1 . Using this approach, we can completely eliminate 
the effect of distance errors in the data being analyzed. 
Indeed, when using the method with known distances 
to stars, we must multiply the left-hand and right-hand 
sides of Eqs. (1) and (2) by 4.74r and r, respectively; 
the unknowns being determined will then be distorted 
by errors in the distances to stars. At present, reliable 
distances to individual stars (with errors < 10%) allow 



us to analyze a solar neighborhood ~100 pc in radius, 
which is not enough for our purposes. 

Let us also consider a Galactocentric cylindrical 
(R, 9, z) coordinate system specified as follows: the z 
axis is directed toward the North Galactic Pole from 
the Galactic center; the azimuthal angle 9 is measured 
from the x axis to the y axis counterclockwise, and R 
is the Galactocentric distance of a star. In the cylindri- 
cal coordinate system, the local deformation and local 
rotation tensor components are 
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provided that the derivatives are taken at point 
{Ro,e o ,z ) = (fl o ,0°,0). 

2 TYCHO-2 

The results of solving the system of equations (1) and 
(2) using Tycho-2 stars of mixed spectral composition 
as a function of magnitude are presented in Fig. 1. We 
used almost all stars of the catalog with the following 
constraint imposed on the absolute value of a star's 

tangential velocity: \fit\ = \J Ha cos S 2 + fig < 300 mas 
yr _1 . The stars were divided into magnitude ranges in 
such a way that each of them contained approximately 
the same number of stars («250 000). In each magni- 
tude range, the random errors of all sought-for param- 
eters are «0.05 mas yr _1 ; the errors in (M+ — M 2 +) 
and (M33 — M+) are twice as large. As can be seen 
from Fig. 1, the parameters that describe the defor- 
mation in the yz plane and the rotation around the 
x axis, i.e., M 2 + and M~ 2 , are almost equal to ze- 
ro. The parameters that describe the deformations in 
the xy and yz planes and the rotations around the z 
and y axes differ significantly from zero. The parame- 
ter M~ is nonzero up to the faintest Tycho-2 stars. 
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The magnitude dependence of M 13 agrees well with 
that found using the TRC catalog, where the mean 
is M- = -0.86 ± 0.11 mas yr" 1 (Bobylev 2004a). 
For Tycho-2 stars fainter than 8 m .5, the mean val- 
ue (the mean often points) is M~ 3 = —0.88 ± 0.10 mas 
yr -1 . This is because the fraction of stars belonging to 
the Local system or local nearby streams is large even 
among the faintest (« 13™) Tycho-2 stars. 

We obtained the following Galactic rotation pa- 
rameters from stars fainter than 8 m .5: M+ = 2.93 ± 
0.21 mas yr -1 and M~ = -2.35±0.10 mas yr" 1 . In this 
case, the Oort constants are A = 4.74M+ = 13.88 ± 
0.98 km s' 1 kpc" 1 and B = 4.74M" = -11.13 ± 0.47 
km s _1 kpc -1 . These Oort constants serve as an indi- 
cator of the consistency of our approach (constraints, 
sky coverage). On the whole, they are in agreement, for 
example, with the results of analyzing the proper mo- 
tions of ACT/Tycho-2 stars (Oiling and Dehnen 2003), 
where only equatorial stars (|6| < 5.73°) were used. 

We use the proper motions of Tycho-2 stars to 
check whether nonuniformities in the latitudinal dis- 
tribution of star s affect the determination of the 
Ogorodnikov-Milne model parameters. For this pur- 
pose, we averaged the stellar proper motions in fields 
of equal area that covered the entire sky. In contrast 
to 48 fields used in the well-known Charlier method 
(Ogorodnikov 1965), we divided the sky into 432 fields. 
The essence of the method is that, despite the differ- 
ence in the number of stars, a unit weight is assigned to 
each field when solving the system of conditional equa- 
tions (1) and (2). As a result, the random errors of the 
main sought-for parameters in each magnitude range 
were ~ 0.15 mas yr -1 ; no marked differences from the 
previous approach were found. The parameter M 2 ~, for 
which the values calculated by this method are indi- 
cated in Fig. 1 by the solid line, constitutes an excep- 
tion. As can be seen from the figure, all of the points 
found have approximately the same displacement of 
w 0.2 mas yr -1 along the coordinate axis. The mean 
Oort constant B is — 12.13±0.63 km s _1 kpc -1 . Appli- 
cation of the method with fields showed that the cor- 
relation coefficients between all of the unknowns being 
determined do not exceed 0.1, except the correlations 
between two unknowns, (M+ — M+) and (M 3 + — M+), 
for which the correlation coefficient is 0.5. Below, we 
use the method with individual stars and will return 
to our comparison of the results obtained by the two 
methods. 

3 UCAC2 

3.1 The First Approximation 

The UCAC2 astrometric catalog is currently the only 
mass catalog that contains highly accurate proper mo- 
tions of faint stars (fainter than 12 m .5) in much of the 
celestial sphere. There is full sky coverage in the decli- 
nation zone —90° < S < 40° and partial sky coverage 



up to S — 52°. The accuracy of the proper motions of 
faint UCAC2 stars given by the authors of the cata- 
log lies within the range 4 to 7 mas yr -1 . The catalog 
extends the ICRS/Tycho-2 system to faint stars. 

Since the catalog contains no observations at de- 
clinations S > +60°, we do not consider the southern- 
most stars at declinations 6 < —60° in order that the 
sky coverage be symmetric. 

The results of solving the system of equations (1) 
and (2) using UCAC2 stars of mixed spectral composi- 
tion are presented in Fig. 2 as a function of magnitude. 
We used a constraint on the absolute of the tangen- 
tial stellar velocity, \nt\ < 300 mas yr -1 . There were 
^250 000 stars in each magnitude range. In the range 
12 m .05— 16 m .01, the stars were taken selectively. In the 
ranges of magnitudes brighter than ll m .17 and fainter 
than 16 m .34, we used all stars from the catalog. Thus, 
we consider a total of ~4.4 million UCAC2 stars. As 
can be seen from Fig. 2, the parameters that describe 
the deformations in the xy, yz, zx planes and the rota- 
tion around the z and y axes differ significantly from 
zero. The parameter M~ 3 is nonzero for the faintest 
UCAC2 stars. 

As can be seen from Fig. 2, there is a noticeable 
jump in the magnitude dependence of M~ 3 for stars 
fainter than 13™, which is inconsistent with the as- 
sumption that this parameter monotonically tends (in 
absolute value) to a minimum. This jump can be ex- 
plained by the presence of a magnitude equation in 
the proper motions of UCAC2 stars. Analysis of Eqs. 
(1) and (2) in the equatorial coordinate system (du 
Mont 1977, 1978) leads us to conclude that the rotation 
around the Galactic y axis can be affected by a system- 
atic error only in /i a cos5. Our simulation shows that 
the coefficient of the magnitude equation for the entire 
catalog is 0.3—0.4 mas yr -1 mag -1 for stars fainter 
than 13 m . Eliminating this magnitude equation yields 
M~ « —1 mas yr -1 for 13 m stars. 

On the other hand, we can study the magnitude 
equation of the UCAC2 catalog by an independent 
method, by comparison with the original Pulkovo data. 
As Bobylev et al. (2004) showed, the Pulkovo stellar 
proper motions do not have a noticeable magnitude 
equation. 

3.2 The Magnitude Equation in UCAC2 

To find the magnitude equation in the proper motions 
of UCAC2 stars, we used the new proper motions of 
faint stars calculated using the stellar coordinates from 
the Pul-3 catalog (Khrutskaya et al. 2004) and current- 
ly available astrometric catalogs: M2000 (Rapaport et 
al. 2001), CMC13 (2003) and UCAC2. 

The Pul-3 astrometric catalog was compiled at 
the Pulkovo Observatory by measuring photograph- 
ic plates taken with a normal astrograph as part of 
A. N. Deutsch's plan to determine the absolute prop- 
er motions of stars in fields with galaxies. The catalog 
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contains equatorial coordinates for more than 50 000 
faint stars to 16™. 5 in 146 fields north of S = —5° for 
the mean epoch of the Pulkovo observations (1963.25) 
in the ICRS/Tycho-2 system. The positional accuracy 
of faint stars in the catalog is 80 mas. 

The M2000 (2.3 million stars) and CMC13 (36.2 
million stars) catalogs were compiled using CCD ob- 
servations with the Bordeaux and La Palma auto- 
mated meridian circles. The two catalogs realize the 
ICRS/Tycho-2 system in the declination zones 11° — 
18° and —3° — 30°, respectively, and contain stars to 
16 m — 17 m . The positional accuracy of faint stars in 
these catalogs is, on average, 40—60 mas at epoch 2000. 

Mutual identifications of stars in these catalogs 
made it possible to compile a new version of the cat- 
alog of proper motions for 34 000 stars. We designate 
the new version of the catalog as PUL3SE. Owing to 
the high positional accuracy of the stars in the original 
catalogs and the large epoch difference between Pul-3 
and the catalogs for epoch 2000, the estimated inter- 
nal accuracies of the proper motions (e AlaC os<5 h e Ma ) lie 
within the range 2 to 4 mas yr -1 . 

For each star included in both catalogs, we formed 
the PUL3SE-UCAC2 differences between the proper 
motions in both coordinates. The mean differences for 
all data are 1 and —0.5 mas yr -1 for the proper motions 
in right ascension and declination, respectively. The er- 
rors of a single difference were found to be <r^ a cos 5 = 
6.3 mas yr -1 and a^ s = 5.7 mas yr -1 . For the groups 
of stars formed as a function of magnitude U2Rmag at 
0"\25 steps, the ratios (e£ o cos5 +e^ Q cosS „ C a c2 )/ CT L cos5 

and ( e w + e %ucac2)/°ls are close t0 unit y- This con- 
firms that the formal estimates of the accuracies are re- 
alistic for both the new proper motions and the UCAC2 
proper motions presented in this catalog. 

The dependences of A/z on magnitude U2Rmag are 
presented in Fig. 3. Analysis of these dependences re- 
veals no magnitude equation in the proper motions of 
faint stars in both coordinates for relatively bright stars 
(brighter than 12 m .5). This may be interpreted 
result of the fact that all of the catalogs used to cal- 
culate the proper motions realize the ICRS/Tycho-2 
system. The Tycho-2 catalog is known to be complete 
to 11m and to contain a considerable number of star s 
to 12 m .5. For fainter stars, systematic magnitude de- 
pendent errors in the coordinates show up in different 
catalogs. 

On the whole, we may conclude that the magnitude 
equation for the proper motions of faint stars in dec- 
lination is indistinct (about —0.5 mas yr -1 ) and may 
be disregarded in our kinematic analysis at the stage 
of our studies in question. 

As can be seen from Fig. 3, the magnitude equa- 
tion for the proper motions in right ascension in the 
magnitude range 13 m — 15"\5 lies within the range 1 
to 2 mas yr^ 1 . In this magnitude range, the variation 
in differences can be represented by a linear law with 
a coefficient of —0.60 ± 0.05 mas yr -1 mag -1 . 



The presence of systematic errors in the coordi- 
nates of stars in the catalogs of early epochs can be 
one of the main reasons for the existence of a magni- 
tude equation for faint stars. Studies by Khrutskaya 
and Khovritchev (2003, 2004) showed that the YS4.0 
catalog (Yellow Sky plates) , which was used as the first 
epoch to derive the UCAC2 proper motions, could be 
the the source of the magnitude equation. As was noted 
in the description to the UCAC2 catalog, the Yellow 
Sky plates were remeasured, but the results have not 
yet been published. 

Platais et al. (1998) showed that there is also a 
significant magnitude equation with a coefficient of 
~ 1 mas yr -1 mag -1 in the published proper mo- 
tions of NPM stars (blue plates) /i a cos S. In con- 
trast to our case, the magnitude dependence of the 
Hipparcos— NPM differences is positive in sign. 

3.3 Solution with the Elimination of the 
Magnitude Equation 

We again solved Eqs. (1) and (2) for UCAC2 mag- 
nitudes fainter than 13 m .0. The proper motions of 
stars /i Q cos<5 with declinations S > — 5° were correct- 
ed for the magnitude equation with a coefficient of 
— 0.60±0.05 mas yr -1 mag -1 . The magnitude equation 
for 16 m .5 stars was assumed to be zero. Since there is 
a mixed effect of the first epochs of the various sources 
used to derive the proper motions in the magnitude 
range 11.5 — 13 m , we do not consider these stars. 

The derived kinematic parameters are indicated in 
Fig. 2 by the solid line and open circles. As can be seen 
from the figure, including the magnitude equation af- 
fected the three rotation tensor components. Compari- 
son of Figs. 1 and 2 leads us to conclude that, first, the 
magnitude dependence of M 2 ~ (the rotation around the 
z axis) for UCAC2 stars corrected for the magnitude 
equation is in better agreement with that found from 
Tycho-2 stars and, second, the magnitude dependence 
of M~ 3 (the rotation around the y axis) corrected for 
the magnitude equation is monotonic, as we expect- 
ed. We attribute the small jump for stars fainter than 
16m to the fact that a transition to distant stars oc- 
curs. The magnitude dependence of M~ 2 (the rotation 
around the x axis) corrected for the magnitude equa- 
tion changes sign, but its absolute value is close to zero. 
At the same time, the deformation in the yz plane un- 
dergoes noticeable variations with magnitude with an 
amplitude of w 0.6 mas yr -1 ; for the faintest stars, 
M+ w 0.3 ±0.05 mas yr" 1 . 

Eliminating the magnitude equation for 13 m stars 
changes M~ by 25% and M~ by 100%. 

On the whole, we may conclude that the proper 
motions of UCAC2 stars corrected for the magnitude 
equation in the magnitude range 13 — 15 m became sys- 
tematically closer to the Tycho-2 system. 
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Ta6jlim;a 1: Local rotation tensor components M 13 and M 32 calculated using distant Tycho-2 (solutions 1-3) and UCAC2 
(solutions 4-8) stars based on the simultaneous solution of Eqs. (1) and (2). 



No. 




d w , pc 


B-V 


V 


J-K s 


m 


M 13 , mas yr 1 


1 

M 32 , mas yr 


1 


89089 


637 ± 55 


< 0.2 


10.17 






-0.488 ± 0.046 


-0.247 ± 0.047 


2 


228230 


578 ± 58 


1 - 1.4 


10.32 






-0.455 ± 0.033 


-0.251 ± 0.033 


3 


130938 


873 ± 100 


> 1.4 


10.40 






-0.371 ± 0.035 


-0.266 ± 0.035 


4 


319129 


641 ± 60 






> 0.5 


16.41 


-0.597 ±0.032 


-0.137 ± 0.031 


5 


269675 


643 ± 63 






> 0.5 


16.55 


-0.486 ±0.033 


-0.177 ± 0.033 


6 


192852 


747 ± 76 






< 0.5 


16.41 


-0.487 ±0.030 


-0.014 ± 0.030 


7 


156702 


774 ± 87 






< 0.5 


16.55 


-0.486 ±0.032 


-0.058 ±0.033 


8 


211086 


975 ± 127 






> 0.8 


10.31 


-0.493 ± 0.045 


-0.018 ± 0.043 


Mean 














-0.485 ±0.023 


-0.140 ±0.036 


1 




486 ± 116 


< 0.2 


10.17 






-0.25 ± 0.17 


-0.27 ± 0.17 


2 




470 ± 90 


1 — 1.4 


10.32 






-0.50 ± 0.11 


— 0.17 ± O.ll 


3 




749 ± 196 


> 1.4 


10.40 






-0.38 ±0.13 


-0.28 ±0.13 


4 




604 ± 126 






> 0.5 


16.41 


-0.45 ±0.10 


-0.35 ±0.11 


5 




651 ±138 






> 0.5 


16.55 


-0.42 ±0.10 


-0.39 ±0.10 


6 




744 ± 177 






< 0.5 


16.41 


-0.35 ±0.11 


-0.13 ±0.11 


7 




754 ± 177 






< 0.5 


16.55 


-0.37 ±0.10 


-0.18 ±0.11 


8 




967 ±292 






> 0.8 


10.31 


-0.32 ±0.13 


-0.22 ±0.14 


Mean 














-0.39 ±0.02 


-0.25 ±0.04 



Note: No. is the solution (or sample), N* is the number of stars, solutions 1-8 in the upper and lower parts of the table 
were obtained, respectively, from individual stars and from the same data with a division into Charlier fields. 



3.4 Estimating the Group Distances 

To estimate the distances to stars, we use a sta- 
tistical method (Oiling and Dehnen 2003). As the 
peculiar solar velocity relative to the local stan- 
dard of rest, we take the values from Dehnen and 
Binney (1998): (Uq,V q ,W q ) = (10.00, 5.25, 7.17) ± 
(0.36,0.62,0.38) km s _1 . We calculate the parallaxes 
using the formulas 

4.74 • X Q 4.74 • Z Q 

Uq Wq 

where X Q and Z Q are the stellar group velocity com- 
ponents that we found by solving Eqs. (1) and (2), in 
mas yr -1 . Since the Yq component is distorted appre- 
ciably by the asymmetric drift (Dehnen and Binney 
1998), this component is not used to determine the 
group parallaxes. We find the distance d from the re- 
lation d — 1/tt. In this case, the distance error can be 
estimated from the relation 




we estimate e n for the motion along the z axis as 




A similar relation can be derived for the motion along 
the x axis. Practice shows (Dehnen and Binney 1998; 
Oiling and Dehnen 2003) that the z components of 



the solar velocity are more stable than the x velocity 
components for stars of various Galactic subsystems; 
therefore, the parallaxes ir\y are preferred to irjj. 

First, we apply the method to the Hipparcos cat- 
alog. The division into groups is performed using 
trigonometric parallaxes of the catalog with e^/ir < 1. 
In Fig. 4, our statistical distances to stars are plotted 
against trigonometric distances. We see from the fig- 
ure that there are no significant systematic distortions 
in the statistical distances and that the errors in the 
statistical distances at distances of more than 200 pc 
are smaller than those in the trigonometric distances. 
In particular, we estimated the mean distance for 4216 
Hipparcos stars with negative trigonometric parallaxes 
and obtained the following values: du — 838 ± 244 pc 
and dw = 716 ±211 pc. This confirms our assump- 
tion that such stars are, on average, far from the Sun 
(Bobylev 2004a). 

Next, we use the method to analyze UCAC2 stars. 
The proper motions of northern-sky stars fainter than 
13 m .O were corrected for the magnitude equation 
found. In Fig. 5, the statistical distances to stars cal- 
culated using the parallaxes ttw are plotted against 
the distances found using nu- As can be seen from the 
figure, there is a difference in the distances 1/irw and 
I/7177 at ~600 pc that gradually disappears near 700 pc. 
This difference corresponds to the stars in the magni- 
tude range 12 — 13 m and stems from the fact that the 
magnitude equation (which was not compensated for 
in this range) affects the determination of the X Q ve- 
locity component. We see that stars fainter than 16 m .O 
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Ta6jiHiia 2: Cross-correlation coefficients derived when simultaneously solving the system of conditional equations (1) 
and (2) using individual stars (above the diagonal) and by the method with Charlier fields (below the diagonal) 



Parameter 




1 


2 


3 


4 


5 


/? 




7 


8 


9 


10 


11 


A© 


1 


1.0 


0.01 


-0.01 


0.01 


0.0 


-0.05 


0.02 


0.15 


-0.04 


-0.01 


0.0 


Ye 


2 


0.0 


1.0 


-0.01 


0.15 


0.05 


0.0 


0.20 


0.03 


-0.01 


-0.04 


0.0 


Zq 


3 


(J. (J 


U.U 


l.U 


-U.U1 


U.U 


-U.1U 


U.U 


U.U 


n no 

U.U8 


n no 


-U.U5 


M + 

21 


A 


u.u 


u.u 


u.u 


l.U 


n no 
-u.uz 


U.U1 


n o^ 

U.UO 


u.u 


u.u 


n ni 

U.Ul 


(1 
u.u 


^32 


5 


0.0 


0.0 


0.0 


0.0 


1.0 


0.0 


-0.02 


0.02 


0.0 


-0.62 


0.03 


M 13 


6 


0.0 


0.0 


0.0 


0.0 


0.0 


1.0 


-0.02 


0.01 


0.61 


0.0 


0.0 




7 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


1.0 


0.01 


-0.02 


0.02 


0.0 


^l + l-22 


8 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


1.0 


0.0 


-0.02 


0.39 


M+ 


9 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


1.0 


0.0 


0.03 


^2 + 3 


10 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


1.0 


-0.01 


1V1 33-22 


11 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.50 


0.0 


0.0 


1.0 



Note: The coefficients are given for sample No. 3 from Table 1. 



are, on average, farther than 650 — 700 pc, with the er- 
ror in the group distances being no larger than 14%. 
The distance found shows that these stars are locat- 
ed, on average, outside the local (Orion) spiral arm or 
the Local system of stars. This gives reason to consid- 
er the rotation around the Galactic y axis found using 
the faintest stars (Fig. 2) as a residual rotation of the 
ICRS/Tycho-2 system relative to the extragalactic co- 
ordinate system. 

3.5 Estimating the Residual Rotation of the 
ICRS/Tycho-2 System 

Results of the simultaneous solution of Eqs.(l) 
and (2). Let us estimate the limit to which M~ 3 tends 
as the distance to stars increases. We proceed from the 
fact that, on the whole, distant stars of various spec- 
tral types belonging to various Galactic components 
may not have any rotation around the y axis. At the 
same time, we consider the presence of a significant 
rotation as a residual rotation of the ICRS/Tycho-2 
system relative to the extragalactic coordinate system. 

To solve the formulated problem, we use almost all 
of the Tycho-2 stars fainter than 8 m .5 and all of the 
UCAC2 stars fainter than 16 m .34 (about one million 
stars). We gave preference to the Tycho-2 catalog in 
the range of bright stars because of its full sky coverage. 
In the UCAC2 catalog, we selected the range with a 
small magnitude equation that we disregard. 

To extend the distance scale where possible, we re- 
ly on a well-known effect related to the fact that stars 
of different brightnesses and colors are, on average, lo- 
cated at different distances. As a result, we divided the 
Tycho-2 stars into two parts in magnitude with the 
boundary V = 11"\25, into bright and faint ones, and 
then divided each part into seven parts in B — V color 
in the following ranges: <0.2, 0.2-0.4, 0.4-0.6, 0.6-0.8, 
0.8-1, 1-1.4, >1.4. 

We use only those stars from the UCAC2 catalog 



for which the photometric magnitudes from the 2MASS 
survey are given, J,H,K S . In this case, it is impor- 
tant to ensure that the stars cover the entire sky. As 
a result, we divided the stars into four parts in both 
magnitude and J — K s color (with the m .5 boundary). 

The results obtained are presented in Fig. 6 and 
Table 1. As we see from Fig. 6, the rotation, M~, and 
deformation, M+, tensor components undergo notice- 
able and coherent variations with increasing distance. 
Analysis of the plots leads us to conclude that there 
is an actual rotation of the stars confirmed by defor- 
mations up to a distance of ~600 pc. This effect stems 
from the fact that the stars belong to the Local sys- 
tem. This boundary is most distinct in Fig. 6b, where 
M+ decreases sharply to zero at distances larger than 
600 pc. 

Table 1 lists the following parameters: the solution 
(or sample) number, the number of stars, the distance 
dw , the B — V and V ranges for Tycho-2 stars, J — K s 
and the mean magnitude (U2R) for UCAC2 stars; the 
last two columns give the derived components of the 
vector of rotation around the y and x axes. Solution 
No. 1 - 3 pertain to bright (8 m .5 - ll m 25) Tycho-2 
stars; solution No. 4 — 8 were obtained from UCAC2 
stars. The solutions obtained using individual stars and 
the same data, but with a division into Charlier fields, 
are given in the upper and lower parts of the table, 
respectively. The maximum number of fields is 432 for 
Tycho-2 and 367 for UCAC2. For each method, the ta- 
ble gives the corresponding weighted mean parameters 
M~ and M+. 

Table 2 gives the correlation coefficients obtained 
when simultaneously solving the system of conditional 
equations (1) and (2) for sample No. 3. As we see, 
applying the method with fields reduces significantly 
the correlation between the unknowns (M~)- (M+) 
and (M 3 ~)- (M 2 +) being determined. 

As can be seen from the upper part of Table 1, M~ 3 
obtained from bright Tycho-2 stars with B — V > 1.4 
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(solution No. 3) is the minimum one (in absolute val- 
ue). Together with distant M-type giants, sample No. 3 
also includes nearby main-sequence dwarfs. Giants can 
be separated more reliably from dwarfs using 2MASS 
photometry (see, e.g., Majewski et al. 2003; Babusiaux 
and Gilmore 2005). We decided to check this result 
using bright (U2Rmag: 8 m .0 — ll m .17) stars from the 
UCAC2 catalog, where the magnitude equation is neg- 
ligible (solution No. 8). For this purpose, we selected 
211 086 stars according to the criterion J — K s > 0.8; 
the stars lie in the K s magnitude range 3 — 7 m (the 
mean is 6 m .75).We chose the boundary J — K s > 0.8 
empirically so as to have a representative sample. Most 
importantly, s=s 50% of the stars lie outside the zone of 
the Galactic plane, i.e., have \b\ > 10°. We have the 
opposite picture for other distance indicators, distant 

and B stars (solution No. 1), since almost all of them 
are distributed in the Galactic equatorial plane, which, 
in particular, reduces the reliability of the M~ 3 esti- 
mate. For solution No. 1, the correlation coefficients 
between the unknowns (M-)-(M+) and (M 32 )-(M+) 
being determined are « 0.85. Comparison of our esti- 
mates of the distances djj and dw for bright UCAC2 
stars with the photometric distances from Babusiaux 
and Gilmore (2005) for the (K s ) — (J — K s ) plane leads 
us to conclude that the estimate of d w 1 kpc for our 
sample No. 8 of stars is quite realistic. 

Applying the method with fields showed that solu- 
tion No. 1 is least reliable. This is because the number 
of stars outside the zone of the Galactic plane (b > 10°) 
is 8 — 10; in all the remaining solutions, the number of 
stars in such fields is an order of magnitude larger. So- 
lutions No. 3 and No. 8 are most reliable. In the strict 
sense, they are not independent, since the Hipparcos 
and Tycho-2 catalogs were used (Zacharias et al. 2004) 
as the first epoch to determine the proper motions of 
bright (to U2R w 12 m .5) UCAC2 stars. 

The mean Oort constants were found from all of the 
UCAC2 stars fainter than 16 m .34 using the method 
of fields to be A = 12.10 ± 0.58 km s" 1 kpc" 1 and 
B = -10.63 ± 0.46 km s" 1 kpc" 1 . 

The estimate of M~ = -0.37 ± 0.04 mas yr" 1 ob- 
tained from solution No. 3 in the upper part of Table 

1 is in good agreement with M~ = —0.38 ± 0.13 mas 
yr -1 calculated by the method of fields (the lower part 
of Table 1) and with the mean M~ = -0.39 ±0.02 mas 
yr -1 calculated using all of the data in the lower part 
of Table 1. 

The mean M~ 3 calculated using the data in the low- 
er part of Table 1 is fairly stable, depending on the 
constraint imposed on \nt\- Above, we used the cri- 
terion \nt\ < 300 mas yr -1 (in two iterations, with 
analysis of the residuals using the 3a criterion), which 
allows nearby stars with high velocities to be reject- 
ed. We calculated the lower part of Table 1 for vari- 
ous constraints on \nt\i from 300 to 25 mas yr -1 . The 
constraint |/it < 25 mas yr -1 implies that at a mean 
stellar distance of 700 pc, the maximum stellar space 



velocity is 110 km s _1 , i.e., only disk stars fall into the 
sample. The means (sample No. 1 was not considered 
because of the small number of stars at high latitudes 
remaining under these constraints) are 

M~ = -0.33 ± 0.02 mas yr" 1 , (4) 
M+ = +0.05 ± 0.06 mas yr" 1 , 
M~ = -0.26 ± 0.04 mas yr" 1 , 
M+ = +0.03 ± 0.04 mas yr' 1 , 
< 25 mas yr -1 . 

Results of the solution of only Eq. (2). Tables 
3 and 4 give the kinematic parameters that were deter- 
mined by solving only Eq. (2). We used two approach- 
es: the solution based on individual stars (the upper 
parts of Tables 3 and 4). For the Tycho-2 catalog (Ta- 
ble 3), we used the samples of stars corresponding to 
solutions No. 2 and 3, whose parameters are listed in 
Table. 1. We imposed the constraint \p,t\ < 100 mas 
yr^ 1 on the absolute value of the tangential stellar ve- 
locity. Using this constraint reduced considerably the 
random errors of the unknowns in the method with 
fields and affected noticeably sample no. 8: the mean 
distance to the stars increased to dw ~ 1100 pc. 

For the UCAC2 catalog (Table 4), we used the 
faintest stars (the stars whose parameters are listed 
in Table 1 for solutions nos. 4 and 6 as well as for so- 
lutions nos. 5 and 7 were combined into two groups); 
about 450 000 stars were used in each magnitude range. 

The results of solving Eq. (2) obtained by the two 
methods show that the Oort constant A and the lo- 
cal rotation tensor components are determined reliably 
and confirm the results that we obtained by simulta- 
neously solving the system of conditional equations (1) 
and (2). 

Table 5 gives the correlation coefficients derived 
when solving Eq. (2) for sample no. 3. As we see from 
Table 5, applying the method with fields reduces the 
correlation between the unknowns (M~)-(M+) and 
(M 3 ~)-(M+) being determined by an order of magni- 
tude. Comparison of Tables 5 and 2 leads us to con- 
clude that the correlation between the above unknowns 
is lower when Eqs. (1) and (2) are solved simultaneous- 
ly; therefore, the simultaneous solution is preferred. 

In our model, we did not abandon the unknowns 
(M+ - M+) and (M+ - M+). The kinematic mean- 
ing of the first difference is that 0.5(M+ - M+) = C, 
where C is the Oort constant. We can determine the 
vertex deviation l xy in the xy plane from the relation 
tan 2/^ = —C/A by analyzing the Oort constants A 
and C. Our results show that the vertex deviation for 
the most distant stars under consideration (sample no. 
8) differs significantly from zero; it is l xy = 7± 1° and, 
for sample no. 3, l xy = 6 ± 2°. 
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Ta6jiHD,a 3: Kinematic parameters calculated from distant Tycho-2 stars using only Eq. (2). 



V 


dw, PC 


A, km s 1 kpc 1 


M 13 , mas yr 1 


M-jtj, mas yr 1 


M 32 , mas yr" 1 


M+, mas yr" 1 


10.32 
10.40 


572 ± 55 
856 ± 94 


13.8 ±0.4 

13.9 ±0.5 


-0.30 ±0.05 
-0.27 ±0.06 


-0.09 ±0.06 
-0.15 ±0.07 


-0.20 ±0.05 
-0.20 ±0.06 


-0.05 ±0.06 
-0.07 ±0.07 


10.32 
10.40 


484 ± 82 
714 ± 165 


12.5 ±1.0 
11.1 ± 1.2 


-0.33 ±0.09 
-0.29 ±0.11 


-0.36 ±0.17 
-0.37 ± 0.21 


-0.19 ±0.09 
-0.22 ±0.11 


±0.16 ±0.17 
+0.04 ± 0.21 



Note: The solutions obtained from individual stars and from the same data with a division into Charlier fields are given 
in the upper and lower parts of the table, respectively. 



Ta6jiHD,a 4: Kinematic parameters calculated from distant UCAC2 stars using only Eq. (2). 



m 


d w , pc 


A, km s 1 kpc 1 


M 13 , mas yr 1 


M-jtj, mas yr 1 


M 32 , mas yr" 1 


M^, mas yr" 1 


16.41 
16.55 
10.31 


731 ± 65 
727 ± 68 
1122 ± 128 


12.3 ±0.4 

13.4 ±0.4 
13.9 ±0.4 


-0.41 ±0.04 
-0.42 ±0.04 
-0.17 ±0.05 


0.05 ±0.04 
0.06 ±0.05 
0.26 ±0.06 


-0.14 ±0.04 
-0.20 ±0.04 
-0.17 ±0.06 


-0.28 ±0.11 
-0.31 ±0.06 
±0.19 ±0.06 


16.41 
16.55 
10.31 


701 ± 142 
719 ± 149 
1184 ± 274 


10.9 ± 1.4 
11.2 ± 1.4 
12.2 ± 1.0 


-0.32 ±0.10 
-0.29 ±0.11 
-0.27 ±0.08 


0.26 ±0.15 
0.35 ±0.15 
0.05 ±0.11 


-0.16 ±0.11 
-0.26 ±0.11 
-0.12 ±0.09 


-0.12 ±0.15 
-0.15 ±0.15 
±0.18 ±0.11 



Note: The solutions obtained from individual stars and from the same data with a division into Charlier fields are given 
in the upper and lower parts of the table, respectively. 



4 DISCUSSION 

On the whole, we may conclude that the estimate 
M~ 3 = —0.37 ± 0.04 mas yr -1 can be interpreted as 
a residual rotation of the ICRS/Tycho-2 system rel- 
ative to the inertial reference frame. This value is in 
good agreement with M~ = —0.36 ± 0.09 mas yr" 1 
that we obtained by analyzing distant Hipparcos stars 
(Bobylev 2004a) using trigonometric parallaxes to an- 
alyze the residual rotation of the ICRS/HIPPARCOS 
system as a function of distance. Previously (Bobylev 
2004b), based on an independent astrometric method, 
we found the following equatorial components of the 
vector of residual rotation of the ICRS/HIPPARCOS 
system relative to the inertial reference frame: u> x = 
+0.04 ± 0.15 mas yr" 1 , uj y = +0.18 ± 0.12 mas yr" 1 , 
and us z = —0.35 ± 0.09 mas yr" 1 . Transforming the ui z 
component to the Galactic coordinate system yields 
My = M~ = -0.26 ± 0.07 mas yr" 1 , which is also in 
agreement with the above estimate. 

Inaccurate realization of the ICRS/Hipparcos sys- 
tem, in other words, a residual rotation of the 
ICRS/Hipparcos system relative to ICRF, can be one 
of the main causes of the rotation found (Ma et 
al. 1998). Since a linear model of solid-body rota- 
tion (Kovalevsky et al. 1997) was used to link the 
ICRS/Hippacros system to ICRF, we believe it appro- 
priate to use the linear Ogorodnikov-Milne model to 
solve the formulated problem. 

The best solution of the problem of controlling the 
inertiality of the ICRS/Tycho-2 system may be asso- 
ciated with the use of images for extragalactic objects 



in the UCAC2 catalog. As follows from the description 
to the electronic version of UCAC2, the measured im- 
ages of star like extragalactic objects are present. In the 
catalog, they are not marked in any way and are not 
identified; their total number is unknown. As follows 
from our work, for the random errors in the sought- 
for parameters to be —0.05 mas yr -1 , we must have 
-100000 objects. 

On the other hand, we see from Fig. 6 that the vari- 
ations in the rotation and deformation parameters with 
distance are nonlinear in pattern. This is attributable 
to peculiarities of the structure of the Local system 
of stars and larger-scale Galactic structures (the spi- 
ral arms and the disk warp). The possibilities for a 
more detailed study of these peculiarities will emerge 
soon, when the next version of the UCAC2 catalog, 
providing full coverage of the celestial sphere, will be 
published. In this case, the magnitude equation should 
be completely eliminated from the proper motions of 
its stars. 

Our results have a bearing on the problem of study- 
ing the Galactic disk warp (Miyamoto and Soma 1993; 
Miyamoto et al. 1993b; Miyamoto and Zhu 1998). 
Since O B5 stars concentrated near the Galactic plane 
are mainly analyzed in these papers, the correlation 
coefficients between the parameters with subscript 
3 are close to 0.95. To reduce the correlations, the 
Ogorodnikov-Milne model is simplified by making the 
assumptions (in our notation) M+ = — M~ and M~ 2 = 
M+, suggesting that dVn/dz = dVg/dz = 0. Under 
these assumptions, the above terms disappear from Eq. 
(1). In this case, the key equation is Eq. (2), which al- 
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Ta6jiHiia 5: Cross-correlation coefficients obtained when solving Eq.(2) using individual stars (above the diagonal) and 
by method with Charlier fields (below the diagonal) 



Parameter 




1 


n 

I 




6 


4 


5 





7 




8 


9 


1 A 

10 


v 

A 


1 


1.0 


0.0 


-0.02 


0.01 


-0.01 


-0.06 


0.12 


-0.04 


0.0 


-0.01 


Y Q 


2 


0.0 


1.0 


-0.02 


0.14 


0.05 


0.01 


-0.02 


0.0 


-0.03 


-0.03 








n n 
u.u 


n n 
u.u 


1 n 
± .u 


-u.uz 


u.u± 


-U.UD 


n n 
u.u 


u.u^ 


n no 
-u.uz 


-U.UO 


M+ 


4 


0.0 


0.0 


0.0 


1.0 


0.0 


0.03 


0.01 


0.02 


-0.01 


0.0 




5 


0.0 


0.0 


0.0 


0.0 


1.0 


0.0 


0.01 


0.0 


-0.89 


0.03 


Mf 3 


6 


0.0 


0.0 


0.0 


0.0 


0.0 


1.0 


-0.02 


0.88 


0.0 


-0.01 




7 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


1.0 


-0.04 


-0.01 


0.81 


M+ 


8 


0.0 


0.0 


0.0 


0.0 


0.0 


0.02 


0.0 


1.0 


0.0 


-0.01 


M23 


9 


0.0 


0.0 


0.0 


0.0 


-0.02 


0.0 


0.0 


0.0 


1.0 


-0.02 


-M33-22 


10 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.82 


0.0 


0.0 


1.0 



Note: The coefficients are given for sample No. 3 from Table 1. 



lows us to determine the quantities 2M+ = — 2M 13 = 
dV z /dR and 2M~ = 2M+ = (l/R)(dV z /d6). The oth- 
er model of these authors is constructed by assuming 
that M" = M+ = 0, M 3 " = M+. In this case, only 
one quantity, 2M 3 ~, is determined from the solution 
of the system of conditional equations (1) and (2). As 
we showed, using the method with fields, which is al- 
most free from any correlations between the unknowns 
(M-)-(M+) and (M-)-(M+), allows the problem to 
be considered without any preliminary assumptions. 
The parameters in Tables 3 and 4 show that M+, M 3 ~ 
and M 2 + have statistically significant differences, de- 
pending on the method of analysis and the sample. 

Using Galactic stars to solve our problem has a 
number of peculiarities related to the structure of the 
Galaxy. Let us consider some of them. 

(1) The yz plane, rotation around the x axis. Chiba 
and Beers (2000) showed that old (metal-poor) thick- 
disk stars lag well behind the circular Galactic rotation 
velocity, dVg/d\z\ — — 30 ± 3 km s" 1 kpc" 1 ; for halo 
stars, this gradient is almost twice as large. This quan- 
tity expressed in mas yr" 1 is 30/4.74 = 6.3 mas yr" 1 ; 
it can have an effect if the distribution of sample stars 
is asymmetric in z. As can be seen from Eqs. (3), the 
gradient dVe/dz enters into the parameters M+ and 
M~ 2 . Applying the method with fields and using the 
constraint on \^ t \ (solution (4)) allows this effect to be 
minimized. The presence of an appreciable component 
M 3 ~ (solution (4)) should be considered as a manifes- 
tation of a particular effect in the actual kinematics of 
stars, since we consider M~ as the main component in 
the problem of controlling the inertiality of the ICRS 
system as a residual rotation (Bobylev 2004b). The re- 
lationship between distant giants and kinematics of the 
warped Galactic disk may be such an effect. However, 
a detailed analysis of this problem is beyond the scope 
of this paper. 

(2) The xy plane, rotation around the z axis. 
The linear Ogorodnikov-Milne model we used allows 
a parabolic fit to the Galactic rotation curve (Vg(R)) 



near R Q to be obtained. Analysis of the currently avail- 
able data shows that the Galactic rotation curve at 
the solar distance is flat (Popova and Loktin 2005; 
Avedisova 2005; Zabolotskikh et al. 2002). Previously 
(Bobylev 2004c), we determined the following param- 
eters of the angular velocity (lo(R)) of Galactic rota- 
tion by analyzing distant (d = 1 kpc) O and B stars: 
uj q = -28.0 ± 0.6 km s" 1 kpc" 1 , u' = +4.17 ± 0.14 
km s" 1 kpc" 2 , lu'J = -0.81 ± 0.12 km s" 1 kpc" 3 , 
R =7.1 kpc. Based on these data and assuming that 
R — R = 1 kpc, we can see that the effect of the sec- 
ond derivative |0.5(i? - R ) 2 lo'o | = 0.4 km s _1 kpc" 1 
is not significant, is comparable to the random errors 
in the Oort constant A and B. Thus, disregarding the 
second derivative w" gives a contribution to the am- 
plitude of the waves that we see in the upper plots of 
Figs. 1 and 2 ( M+ and M~) of -0.1 mas yr" 1 , which 
is considerably smaller than the amplitude observed in 
the plots. 

5 CONCLUSIONS 

Our analysis showed that the proper motions of 
UCAC2 stars are distorted by a complex magnitude 
equation. Thus, for example, no noticeable effect of the 
magnitude equation was found in the proper motions 
of stars brighter than w ll m , but this effect is signif- 
icant in the magnitude range 12 — 15 m . We showed 
that the magnitude equation has the greatest effect in 
[i a cos S of northern-sky stars. This is because differ- 
ent sources for the southern and northern skies were 
used as the first epoch in deriving the proper motions 
of UCAC2 stars. We determined the coefficient of the 
magnitude equation for northern sky stars fainter than 
13 m by comparing the proper motions of UCAC2 and 
PUL3SE stars, -0.60 ± 0.05 mas yr" 1 mag" 1 . We as- 
sociate this magnitude equation with the effect of the 
NPM catalog, Yellow Sky plates. The derived magni- 
tude equation mainly affects the determination of the 
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rotation tensor components. Its elimination leads to 
a 25% change in the rotation around the z axis (the 
Oort constant B) and to a 100% change in the rota- 
tion around the y axis for 13 m stars. 

The method of statistical parallaxes was used to 
estimate the distances to stars with random errors no 
larger than 14%. The linear solar velocity relative to 
the local standard of rest, which is well determined for 
the local centroid (d ~ 150 pc), was used as a refer- 
ence. This could affect the distance estimates for dis- 
tant stars. 

The method with fields was shown to reduce the 
correlations between the unknowns (M 1 ~)-(M+) and 
(M 3 ~)-(M+), being determined by an order of magni- 
tude. When Eqs. (1) and (2) are solved simultaneous- 
ly, the correlations between these unknowns are lower 
(than those for the solution of only Eq. (2)); therefore, 
the simultaneous solution is preferred. 

The mean Oort constants estimated from Tycho- 
2 stars (8 m .5 - ll m .25) are A = 13.88 ± 0.98 km 
s" 1 kpc" 1 and B = -12.13 ± 0.63 km s" 1 kpc" 1 . 
The analogous means that we estimated from bright 
(8 m - ll m .17) UCAC2 stars have no statistically sig- 
nificant differences. The mean Oort constants estimat- 
ed from UCAC2 stars fainter than 16 m .34 are A = 
12.10 ± 0.58 km s" 1 kpc" 1 and B = -10.63 ± 0.46 km 
s- 1 kpc- 1 . 

We established that the two components of the 
Ogorodnikov-Milne model that describe the rotation 
around the Galactic y axis and the deformation in 
the xz plane depend on heliocentric distance. For the 
nearest stars (d w 200 pc), these parameters are 
— 1.0 ±0.1 mas yr" 1 and decrease, on average, by a fac- 
tor of 2 at distances of ~600 roc. We associate this with 
kinematic peculiarities of the Local system of stars. For 
more distant stars, there are, on average, no statisti- 
cally significant deformation in the xz plane. 

For distant Tycho-2/UCAC2 stars at heliocentric 
distances of «900 pc, the mean rotation around the 
Galactic y axis is —0.37 ± 0.04 mas yr -1 , which we 
interpret as a residual rotation of the ICRS/Tycho-2 
system relative to the inertial reference frame. 
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Fig. 1. Kinematic parameters inferred from the proper motions of Tycho-2 stars vs. magnitude. 
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Fig. 2. Kinematic parameters inferred from the proper motions of UCAC2 stars vs. magnitude. 
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Fig. 3. PUL3SE-UCAC2 stellar proper motion differences vs. magnitude: (a) A^ Q cos(5, (b) Afig. 
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Fig. 4. Statistical distances to stars vs. trigonometric distances calculated using Hipparcos data. The division 
into groups was performed using trigonometric parallaxes. 
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